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IBSTMCT 


A one-dimensional water quality simulation model which 
incorporates mass transport by longitudinal dispersion 
and advection is presented. The model simulates temperature 
dissolved oxygen, and conservative minerals. Dissolved 
oxygen balance has been coupled with carbonaceous bio- 
chemical oxygen demand, nitrification and atmospheric 
reaeration. The model was used to simulate water quality 
parai.ieters in tlahi river during relatively dry period of 
the month of May. Simulated results matched fairly well 
with the observed values. 



CHAPTER 1 


lETRODUCTIOH 


Rapid urbanization and industrialization has left mankind 
with problems of waste disposal. liq,uid wastes generated 
in a community and industrial liquid wastes generated in 
industries are disposed off in natural water bodies like 
strecims , rivers or sea. As a result of rapid industrialization 
the waste disposal in rivers has increased and has been 
localized near urban and industrial areas, which leads to 
degradation of water in its physical, chemical and biolo- 
gical qualities, Eor many engineering purposes, the 
concentrations of dissolved oxygen (DO) and decomposable 
organic matter in the stream were used as indicators of 
water quality (Streeter and Phelps, 1925). Because of 
complexities of various industrial wastes discharged into 
rivers, complex reactions and inter-relationship of various 
mechanisms including nitrification, oxygen uptake by oxygen 
demanding chemical wastes, etc. need to be coupled in 
assessment of dissolved oxygen balance in water body. 

For the assessment of the extent of pollution in rivers, 
costly and time consuming water qxoality surveys are being 
under-taken and decisions for anticipated conditions in . 
future arc carried out usually by intution and experience 
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gained for similar conditions elsewhere. Engineers have ' 

constructed models to study such complicated systems. 

Modelling and simulation are techniques frequently used 
(and sometimes ahusedl) in today's scientific and 
engineering investigations. 

A model may he thought of as being a representation of a i 

system in a form suitable for demonstrating the way the 
system behaves, while simulation involves subjecting ; 

models to various changes in such a way as to explore the 
possible effects of these changes on the real system. ¥hen ■ 

properly applied modelling and simulation can result in I 

considerable saving in both time and money. | 

i 

t 

The mathematical model is a mathematical representation of 
the major mechanisms in a natural system in such a form 

I 

that a cause and effect relationship can be analjrtically | 

I 

, I 

approximated. A complicated system like river system (for | 

Water Quality Analysis) may require the combination of 
several mathematical models. These models can be transfoimed 
in such a way that output of one model can be used as input to j 

the other and incorporated into a computer program to form [ 

1 

a computerized simulation model. 

1 • 1 Object iye^ . .PJL jfc » 

The present investigation was undertaken to construct a | 

mathematical model for predicting dissolved oxygen concen™ I 

tration and other water quality parameters in Mahi RLver System, | 
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X'j'hich. receives industrial wastes from industries near 
Baroda. These industries are Gujarat Refinary, Gujarat 
Industrial Development Corporation Estate, Gujarat State 
PertiD.izer Corporation, Indian Petro-Chemical ltd., Suhrid 
Geigy (a pharmaceutical firm), Universal Dye Stuff, etc. 

Rapid development in production in these industries has 
caused increased pollution in river Mahi. National Environ- 
mental Engineering Research Institute (NEERI), Nagpur 
carried out Mahi river survey in 1972-74. The work was 
sponsored by Gujarat State Public Health Department. NEEEI 
proposed to construct an open channel to carry waste from 
industries mentioned above to much down-stream in estuarine 
region of Mahi river (Baroda Effluent Channel Project, Ref. 2) 
It was sought to suggest in the present study the feasibility 
and need for constructing the channel. 

On the Mahi river, there is a multi-purpose dam Kadama. 
Another multi-purpose project is under execution on Panam 
river, a major tributary of Mahi river. It is intended to 
evaluate quantitatively the water-quality oriented alter- 
natives such as waste treatment level of combined waste 
from above mentioned industries versus flow augmentation from 
Kadana and Panam reservoirs in the Mahi river system. 
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LITEMTUEB REVIEW 

Environmental considerations have concerned water resource 
planners for many years. In the last 10 years, public 
interest has caused significantly increased efforts towards 
development of numerical tecliniques for analysis of water 
quality conditions in water resources systems. 

2-1 iu Modelling ; 

Finite element model is one of the common methods in modelling 
river system. Finite element model of a stream system 
consists of a series of elements as shown in Fig. 2,1 
corresponding to a discrete stream segment arranged so 
that the output from one element becomes the input to the 
next. The transfer function is determined by perfoiming a 
mass balance of a given water-quality parameter over a time 
interval, At, on a stream segment of cross-sectional 
area, A and of lengths Ax along x-axis. 

2.1,1 Hydrol o gic Balance ; 

Assuming steady-state conditions, a hydrologic balance for 
the control volume (see F ig. 2,1) can be represented by, 


0 


( 1 ) 
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wliere P is the precipitation rate, E is the evaporation rate* 
Q. 1 is the flow entering i"th control volume, . x is flow 
leaviii/- i"th control volume and is the run-off entering 
or lea-ving i-th control volume, 

2.1,2 Mate r ia 1 B alanc e % 

In a river system, material balance of constituent is carried 
out hy balancing mass entering the control clement, and 
mass leaving the control element. Mass transport of consti- 
tuent is accomplished by advection and dispersion. Disper- i 
' i 

sion represents the combined effects of density differences 
and velocity gradients. Assuming only longitudinal disper- j 
Sion to bo predominant a temporal distribution of a consti- | 
tuont within a control volume due to dispersion is given by, 

H = (2) I 

I 

I 

where c = concGntration of constituent under consideration^ ! 

I 

t = some point in time, x = some point along x-axis , 

! 

= longitudinal dispersion coefficient, A = cross-sectional ; 

t 

area of the element, 

I 

Mass transport due to advection will result in reduction of 
the temporal change in concentration of a constituent. Then 
eq^. 2 will become. 
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A = a (ADj^ 1-“) /' 3X - 3{Auc) /ax (5) 

whore u = mean stream velocity. 

Let 'S' be the sink or source strength of a constituent under 
question. Then mass transport equation will be, 

~ 3x^ ± A'S' (4) 

Sq. 4 is a linear, parabolic partial differential equa tion 
describing temporal and spatial variations in constituent 
due to longitudinal dispersion and advection. 

For steady-state conditions. 


3c 

3T ~ 


0 


3 ( ALj^ /3 X - 3 (Auc ) /a X 


+ A' S' 


(5) 


Equation 5 is a sec. order linear differential equation in 
which general boundary conditions are. 


c(0) = 

o(x) = 0^ 

whero = the concentration of particular constituent 
at head water source and C = the concentration of const! 

X 

tuent just downstream of the last element. 
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Solu'Gion techniq.UG for eq.. 5 is described in Appendix A. 
Lo ngitu d i na l Dispe rsion s 

Dispersion is basically a conser-vative transport associated 
with spatially averaged velocity variation, as opposed to 
'diffusion' which is reserved for transport that is associated 
primarily with time-averaged velocity fluctuations. 

Longitudinal dispersion in. rivers is given by eq. 6 (Ref, 12 ) 
D^ = 22.6 n u (6) 

where D]^ = longitudinal dispersion coefficient ft /sec, 
n = Hanning's roughness coefficient tabulated for different 
types of river channels in Table 1. u = Mean velocity, 
ft/sec, and D = moan depth, ft. 

Table 1: Values of Hanning's 'n' roughness 

Coefficient After Henderson (1966). 

H iver Channel s . n 

1, Clean and straight 0,025 - 0.030 

2, Winding with pools and shoals 0.033 ~ 0.040 

3. Very weedy, winding and overgrown 0.075 - 0.150 

1 /6 

4. Clean straight alluvial channels 0,031 d"^' 


(d = size in ft. ) 
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Typical values of dispersion coefficients are given in 
Table 2 (Ref. 12). 


Table 2; Typical Values of Dispersion Coefficients 


System Classification 


Br- ft^/sec. 


1. Plumes and small streams 

2. Large rivers 

3 . Estuaries 


0.05 - 3.0 
5 - 300 

500 - 5000 


2.1.5 Sour ces an d Sinks : 

Pollowing are the descriptions of sources and sinks for 
various non-conservative constituents as described in eq..4. 

(i) Temperature; 

In the thermal behaviour of a water body, it is essential 
to have a quantitative representation of the heat fluxes 
between the xfater surface and the atmosphere. These heat 
fluxes describe sources and sinks for computing variations 
in temperature. A body of water cools by losing heat to 
the atmosphere? conversely, it warms by gaining heat from 
the atmosphere. All bodies of water cool or discharge heat 
to tho atmosphere by back radiation, evaporation and conduc- 
tion, at the same time warming or receiving heat throu^ 
shor'e-wave solar radiation and long-wave atmospheric 
radiation. 
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It is illustrative to represent the net heat flux at the 
water surface as shown in Pig. 2,2. The range in magnitude 
of monthly average values of heat transfer are representa- 
tive of northern latitudes. 


Hg = Short-wave solar radiation (400 to 2800 Btu/ft.^day) 



li = long-wave atmospheric radiation (2400 to 5200 
, Btu/ft? day) 



= long-wave hack radiation (2400 to 3600 
Btu/ft? day) 



= Heat loss due to evaporation (2000 to 
8000 Btu/ft,^ day) 



Conduction heat losses or gains 
(-320 to +400 Btu/ft^ day) 

= Reflected solar radiation (40 to 
200 Btu/ft.^ day) 


H = Atmospheric reflection (70 to 
A 120 Btu/ft.^ day) 



PPig. 2.2; Net Rate at Which Heat Crosses the Air-Water 
Interface. 

The espression that results from the summation of these 
various fluxes is, 


11 

Hsn 


H = 
sn 

Hs 

- 

II 

Id 


"^ar 


( 7 ) 
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■^.’'herG, 

= Net energy flux passing tlie air-water 
interface Btu/ft.^ day 

= Net short-wave solar radiation fluz passing 
through the interface after losses due to 
absorption and scattering in the atmosphere 
and by reflection at the interface, Btu/ft.^ day. 

H = Not long-wavG atmospheric radiation flux 

SLIl 

passing through the interface after reflec- j 

i 

tion, Btu/ft.^ day. I 

= Outgoing long-hack radiation flux Btu/ft.^ day. I 

Hq = Convective energy flux passing back and forth i 

between the interface and the atmosphere, | 

Btu/ft.^ day, and | 

H. = Bnergj’" loss by e-vaporation, Btu/ft,^ day. i 

’ i 

I 

1 

Various heat fluxes and their transfer mechanisms have been I 

j 

well understood and are adequately documented in the lite- | 

nature by Edinger and Gever (1965) as well as in Eef. 12 j 

1 

A stream may be considered completely mixed in its trans- | 

verse section if the heat exchange with the environments^ | 

affects the water body practically instantaneously over its 
entire depth. Wunderlich (1969) found that the internal 

j 

turbulence in many stream was sufficient to assume practically j 
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instantp.ncous heat dis trihution from the surface downward. 
However, as velocities decrease and depths increase, this 
assumption, "becomes less valid. 

Thus, assuming complete mixing, Eq. 4 can "be written with a 
source term for temperature ?,s, 

A = 3 3 . -- 3('Auc)/3x + A'Sg’ /yCp ( 8 ) 

whore Sjj = and D = Hydraulic depth of the stream, 

Y = weight of water, 62,4 Ib/ft^, Cp = Sp. heat of water 
1,0 Btii/lb „ Op. 

The source term ' > which has units of Btu/ft.^ hr, 

accounts for internal heat generation and all heat trans- 
ferred across the system boundaries. That is, heat trans- 
ferred across the air-water interface and heat conducted 
across mud-water interface. In absence of ground water 
flovr, heat is transported across the mud -water interface 
only by molecular conduction which is relatively insignifi- 
cant in comparison to surface heat exchange. 

( i i ) Dissolved Oxygen : 

Source and/or sink term, as the case may be, for dissolved 
oxygen Was defined by Streeter and Phelps' (1925) as. 



K2 (Cg 


c) - I 


(9) 
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whore c = Dissolvod oxygon concentration in mg/l, 

= Soluhility of oxygen in water, mg/l, = atmospheric 
roacration coefficient day”^, = Carbonaceous BOD decay 
rate, day L = Carbonaceous BOD concentration, mg/l, 

Cg = Solubility of oxygen at given temperature is given by 
(Standard Method, 1971), 

Cg = 24.89 - 0.426 T + 0.00373 - 0.0000133 

where, T = Temperature in 


For ruaeration coefficients, six methods are available. 


vO.969 

(1) Churchill et al (1962), IL^ = 5.051 (10) 

H * 


(2) O'Conner and Dobbins (1958), Kp = 3.951 


0.5 


K 


(5) Owens et al. (1964), 

(4) Langbien and Durum (1967),E2 


V' 


0.67 


2 - jjl.85 


3,133 ( 15 ) 

0.5n tL 


( 12 ) 


(5) Tliackston and Kronkel (1966), Kp= 24.95 (l+F * ) g— (14) 


where, K2 = Reaeration coefficient, V = Stream velocity in 
m/soc,, H = Hydraulic depth in meters, F = Froude number = 
V/(Hg)*^*^, = Shear velocity, m/sec. = (H.S.g)*^*^, 

S = Slope of water surface, m/m, g = acceleration due to 
gravity, m^/sec. 
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VJhen industrial wastes are discharged into river, many 
moclianisms occur in the water “body, which are responsible 
for changes in the dissolved oxygen concentration. For 
example, ammonia discharged into river is oxidised to 
nitrites and which are in— turn oxidised to nitrates. During 
this process of nitrification there is considerable amount 
of oxygen demand. 

Several workers, Thomann, R.V. (1972), Willis et al (1975), 

I?.of, 17, have developed dissolved oxygen models which describe 
most of the mechanisms (physical, biochemical, chemical and ! 

biological) mathematically, which affect dissolved oxygen. | 

I 

I 

Hydraulic Engineering Centre of US Army Corps. (Ref.\-7 ) has ; 

1 

developed the ecological model for dissolved oxygen ±as I 

given below. i 

O 2 HO 2 “ HDET (DET+S) O 2 DET - I BIO.O 2 R I 

1 

[BIOR + BIOS (;-Cj^ - 1 )( 1 -BXF)] I 

+ 2 'a( 02 I‘.AU - O2II.AR) (15) I 

I 

where c = Concentration of dissolved oxygen, K 2 = atmosphe- 
ric reaeration coefficient, C„ = concentration of dissolved 
oxygen at saturation, Kn = Rato of BOD removal by oxygen uptake,! 

-L ' 

L = Concentration of ultimate BOD, = Ammonia decay rate, | 

I 

^ i 

!' 
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= Ammonia concentration as nitrogen, 02^11^ = Stoichio- 
metric ecLui valence between oxygen and ammonia, KNO 2 = nitrite 
decay rate, NO 2 = nitrite concentration as nitrogen, 

02K’02 = Stoichiometric equivalence between oxygen and 
nitrite, KDET = Detritus decay rate, DET = Detritus concen- 
tration, S = Concentration equivalent of organic sediment, 
O 2 DEI' = Stoichiometric equivalence between oxygen and 
detritus decay, BIO = Biota concentration excluding algae, 

O 2 R = Stoichiometric qquivalcnco between oxygen and biomass 
respiration, BIOR = Biota respiration rate, BIOG = Biota 
growth rate, BIEPP = Biota digestive efficiency, EXP = Parti-- 
culato fraction of total excrement, A = Algei concentration 
(i.e., phytoplankton and benthic algae), O 2 P = Oxygenation 
factor for algal photossmthesis, AG = Algal growth rate, 

AR = Algal respiration rate. 

( iii ) Bioche mi cal O. xyge n Bemand ; 

Sink equation for BOD in described as, 

a = - L C16) 

( 1"^) Ammonia Nitrogen ; 

+ KDET (DET+S) DE + 

[BIOR + BIOG - 1)(1- EXE)] 

- j; A.AP (AG.Pra - AR) 


[BIO.BIOE 



( 17 ) 
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where JBIOIT = Nitrogen fraction of biota, AN = Nitrogen 
fraction of algae, PT-IN = Ammonia fraction of available 
nitrogen. Other symbols are same as explained in Eq. 15 . 

( V ) Nit r^t e Nitroge n : 

dNOp 

= KNH^.NH^ - ICNO 2 .NO 2 (18) 

where NO 2 = Nitrite nitrogen concentration 
( vi ) Nit rate Nit roge n i 
dNOv 

-g-rp = KNO 2 .NO 2 - I A. AN. AG (1 - ENN) (19) 

Symbols aro same as explained above. 

2 2 Applicatio n of S i mulat ion Models : 

Most of the individuals or groups developing water quality 
simulation models during last few years have recognised 
their efforts as being only a first step in an evolutionary 
development of their respective models. However, many 
worker's have put in their efforts in application of these 
models for real systems. 

Rutherford and O'Sullivan (1974) have made attempts to 
predict water quality in tke Tarawera river using the 
Streeter-Phelps * and various other first order models 
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but none of these could successfully predict the observed 
concentrations of dissolved oxygen. 

The BoX“ Jenkins method, a time based technique for time 
series analysis is successfully used by Huck and Parquhar (1974) 
to model chloride and dissolved oxygen data for St. Clair 
Eiver near Coruna, Ontario, This is thought to be the first 
application of the method to water quality data. The tech- 
nique is demonstrated to be superior in this situation to 
either a frequency-based approach or a deterministic causa- 
tive model. The description of the model building process 
includes the identification, estimation and diagnostic 
chocking stages. 

A ona-dimensional model is presented by Aiba and Ohtake (1976) 
to simulate the spatial and temporal concentrations of 10^-P 
in shallow and polluted river. The model incorporates, 
other than the convection and dispersion, various physico- 
chemical and biochemical reactions of phosphorous sinks and 
sources. With reference to field data on Tama-gava river, 
which penetrates through metropolitan area of Tokyo, the 
model is confirmed to represent the concentration of PO^-P 
fairly well in the mid-region except for the mountaneous 
origin and its estuary. 
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V/illis , R. ,ct.a,l. , ( 1975 ) have simulated water quality in 
Tiuclcee River System in northern California and Nevada. 

They have incorporated tho interaction of the dissolved 
oxygen resources of the water system with the nitrogen 
cycle, ^chlorophyll a and carbonaceous biochemical oxygen 
demand. Also, conservative substances and phosphorous 
uptake by algae and conforms are simulated by their model. 
Mass transport of constituents is accomplished by advec- 
tion only. 

Knowles and Nakeford (1978) have developed a comprehensive 
river water quality model, deterministic in type, in which 
eleven processes that were represented include nitrification, 
denitrification, photosynthesis, BOD decay, reacration and 
others. This deterministic model has been used to predict 
river qu/rlities in Blaclrwator System (part of the River 
Thames Catchment). It was tested under relatively dry 
conditions of 1973 and more normal conditions of 1974. 





SCOPE OP IViE PRESEIE STUDY 


It ic apparent from th.c litoraturo that much insight has 
boon clovolopod in the field of simulation of water quality- 
in rivers. In India, many of tho important rivers are 
being surveyed for their water quality and enough infor-'ia- 
tion will bo available soon to implement such simulation 
raodols on Indian rivers. 

Tho scope of present study is to adapt to Indian conditions, 
tho avro'.lablc know-how in tho field of simulation of water 
quality. Sufficient data arc available for Mahi river system 
in G-ucarat. An attempt has been made in this study to simu- 
late iiroortant paramoters of -^fator quali-tj!- in Hahi river 
by using’ a suij;'.blo model. 

The interrelation and interdependence of various mechanisms, 
which affect water quality, in aquatic environment are 
numerous. Mathematically it is possible to incorporate 
them in tho simulation models as sources and sinks. However, 
it is difficult to quantify these sources and sinks in a 
real system. Hence;, sources and sinks for various water 
quality constituents, which represent such mechanisms, have 
not boon considered. A few of these sources and sinks are 


as follows. 
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1. Adaorption of vr.rioufj conservative and non-consorvative 
substances on organic matter;, jurbidity- -causing substances 

etc . 

2. Denitrification, phot onjnithcsis, BOD removal by 
sedimentation, benthic oxygen demand , effect of biological 
life on dissolved oxygen balance, d '..pendonce of solubility 
of oxygon on turbidity and chlorides etc. 

Dissolved oxygon balance will bo influenced directly or 
indir;.‘Ctly by whole of aquatic environment . It is difficuD.t 
and beg^'ond the scope of this study to take into account 
the entiru aquatic ecosystem, while modelling such complex 
syst on. 

These are sove.ral alternative strategics for control of p 
pollution of rivers. Some of tliGra arc; 

(l) Augracnta.tion of flovr from stora.ge reservoirs, (2) Treat- 
ment for waste loads before disposing the same into the river, 
(3) AD.location of discharge points for waste loads, (4) In- 
stream rinichanic'M roaeration etc. 

From asiong the above the folJ.owing alternatives wore selected 
for this study: 

(1) Plo^''' augmentation from hadana and Panam resorviors, 

(2) Treatment levels for the wastes being disposed in Mini 



rivor and (3) Diverting the uaete waters through a 
clianno?. to dispose it in -'die estuarine region of 



CHAPTER 4 


METHODOIO&T OE STUDY 


The present study was carried out in three distinct phases 
(i) Collection of data, (ii) Verification of model, (iii) Use 
of model as decision making tool. Before discussing these 
phases, a description of Ilahi river system is given below. 

The Uahi river originates from the Vindhya hills near 
about the village of Sardarpur in Dhar district of Madhya 
Pradesh. It flows through the State of Madhya Pradesh and 
&uja,rat over a distance of about 600 km. before it joins the 
Gulf of Capibay. Some of the principal tributaries are 
Som, Anar, Panam, Earad,. Heshri, Mini and Bhadar rivers. The 
section of Mahi river which is under present study is shown 
in jj'ig. 4.1 and schematic of model is shown in Pig. 4.2. 

The Hahi river is the main source of water supply for the 
industrial zone of Baroda city and a supplementary source 
for the Baroda Municipal Corporation. On an average 25 to 
30 mgd of water is derived from radial wells sunk below the 
bed of the river. 

The industrial waste is disposed of in Mahi river through 
its tributary, Mini- Mini river originates at a distance of 





Source of Mahi river 
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Explanation 

Headwater source 

Waste or tributary inflow 

Withdrawal 

O Sampling site numbers 
corresponds to those on figs 
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^ Stream reach numbers 
7 and poiffts of stream 
^ segmentation used in mDdeling 
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TRlBUTARY-3 

(Mini river which collects 
wastes from several 
Industries near Baroda) 

SITE 6 
SITE 7 



70 
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(E) panam river 


TRIBUTARY-1 


TRIBUTARY-3 


- WITHDRAWAL 
FOR BARODA 
W/S 


Fig.4-2 Schematic of model system 
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30 IcQi. north east of Saroda. Its course is almost parallel 
to tint of ilahi. 

The total of 6 mgd waste water is discharged into Mini 
rive'c from various industries. (Ref. 2). 

4 1 JPX ‘ 

Water quality, hydrologic data and local meteriological 
data were collected from various government organisations. 

4.1.1 Hjd_r^jO£;i c^ J^S'jta ; 

Average stream velocity and average depth at three guaging 
stations as shown in Fig. 4.1 were available for three 
different discharge conditions. These data were supplied by 
rivo:;' guo.ging division of irrigation department, Government 
of Gujarat and Material Testing Division of Gujarat 
Engineering Research Institute, Baroda. Monthly average of 
flow rates in Mahi river at Kadana village just downstream 
of luidana dam and monthly average flow rates in Panam river 
for 1977 were also supplied by Irrigabion department. 

The stage-discharge rating curves are prepared from the 
hydrolic data at three guaging stations. These curves 
(Fig. 4.3) are used to determine other hydraulic charac- 
teristics of the stream. Stage-discharge rating curves at 
Sevalia are used for first three reaches, that at Kanora 
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were used for next three ree.cbes and same at Vasad were 
used for last three reaches. 

4,1.2 Wa t^__ Syi'SJ- Z. > 

V/ater quality survey of tlahi river is being carried out by 
Material Testing Division of Gujarat Engineering Research 
Institute (GBRI), Baroda, Water quality data at eight 
sampling points during 1978, as well as flow rates in 
tributaries were j-aade available from the GERI, Baroda. 

Data regarding qirality of wastes were compiled from Baroda 
Effluent Channel Project Report. (Ref. 2 ). V/ater quality 
and flow rate data for wastes and tributaries have been 
presented in Table 3. 

4*1.3 I' p cal_ Glimato lo gic al D ata ; 

Dry bulb tempera tui’e, wet bulb tempemture, atmospheric 
pressure and wind speed, data were available from MeteriolO' 
gica,l Observatory in Faculty of Science, M.S. University 
of Baroda. Five observations per day during May 1978 of 
above mentioned climatological data are used to route 
temperature in Mahi river. 



Taljle 3; Vater qiiality of tributary and ifaste inflows. 
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4 . 2 yerificat ion of Data ; 

rrcliminary computation was carried out using Streeter 
and Phcli^:!' model for dissolved oxygen balance. It is found 
to be inadequate for lower reaches on the downstream of 
confluence of Ifehi and mini river. 

Honco, model developed by HEC as shown in literature is 
used whic ' includes nitrification also. This model is 
compared with the dissolved oxygen concentrations- observed, 
during i-iay 1978. The model for conservative constituent 
is also tested by comparing 'observed concentrations of 
the constituent. 


4*3 U? del • 

It io conejudod in the second phase of study ths.t model is 
well compatible with the real systom. Hence, three state 
variables, viz, flow augmentation, treatment efficiency of 
mini river waste and reduction in waste discharge from mini 
river arc varied to know their impact on the real system. 

Results of both the phases are presented and discussed 
in the next chapter. 



Chapter 5 


K3SUITS AlTD DISCUSSIONS 

As closcrihod oarlior, 1978, sutnmer data were used to verify 
the nodol, Si“iulatod profiles of dissolved oxygen using 
various atmospheric reaoration coofficient formula© are 
shown in fig. 5.1 through Dig. 5.5. All of those profiles 
have siiAilar characteristics upto river milo 70. Drom river 
milo 70 to river mile 90 roaeration formula given by 
ChuD'Chill pt.g^__. (1962) seem to be woll comparable with 
tho observed concentrations of dissolved oxygen. Hence, 
oquabion for K 2 by Churchill et,.al.. (1962) is used in 

further analysis by using tho model as a decision making 
tool. 

Sudd on increase in dissolved oxygon at river mile 5 is 
bocauoQ of Panam river, which moots Mahi river at this 
point with very high flow rate and high dissolved oxygen 
concentration. 

Similar rise at river mile 70 is observed. At this point 
water in withdrawn for Baroda water supply. Due to this 
reason, there is a significant reduction in depths. Because 
of this decrease, the reaeration coefficient is increased 
and honce, dissolved oxygen ooncentration rises at this point. 




paApSSIQ 


River miles 







l/6ui‘ya6Axo paA]OSSiQ 


Fig.5‘4 DO Simulation results 



|/Buu*uaBAxo paAiossiQ 


Rivers miies 

Fig,5*5 DO Simulation results 
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Pigs. 5.6 and 5,7 show Simula bed curves of conservativo 
minerals and trace olomerLts (Arsenic and Cyanido). All the 
computed profiles have much lower values than oToserved 
concentrations. This discrepancy may bo duo to following 
reasons. First, huuian errors encountered during collection 
and analysis of waste samples at disposal points. Secondly, 
in the model formulation, source and sink term is absent 
for conservative minerals, because it is difficult to 
assess 'the extent and location of removal or introduction 
of such conservative minerals in the river system. The 
discrepancies in the results of these minerals are also 
indicates there could be an unidentified source of pollu- 
tion in bliG system, which have been overlooked due to 
ignorance. 

Oynide is biodegradable at a very slow rate. Hence, assim- 
ing it to be a conservative material, concentration of 
cynide in river was simulated. The descrepancies in results 
of this constituent can be dun to this assumption. 

Phenol was also considered as conservative constituent. But 
it can be adsorbed on clay minerals or organic matter which 
then removed from the system by sedimentation. 


Arsenic concentrations wore also simulated assuming it to 
be conservative. Sinks like uptake of arsenic by fishes, 
which are later removed from the system, have been overlooked. 
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Thus, most of the discrepancies in the verification of 
model can bo contributed to the ignorance of many sinks 
or sources, which are assumed norL“existcnt in the model 
formula 'cion. Also, availability of only limited number 
of observations at limited number of sampling points, 
errors involved in collection and analysis of samples 
are oqual contributors of the total errors involved. 

As explained earlier, the model uses stage-discharge rating 
curves for computing velocity and depth (fig. 4.3) in each 
olomento. Those rating curves are constructed from a very 
limited number of observations and are used over long dis~ 
tancos for computing hydraulics of river. Now, these 
variances in velocities and depths can significantly alter 
roaeration contribution to the dissolved oxygen concentra- 
tion through dependence of the li 2 reaeration coefficient on 
velocity. 

Various heat fluxes, that are transferred from water to air 
and visa versa are shoi-m in Figs. 5.8 and 5.9, The simulated 
curve for temperature also lies, as shown in Fig. 5.9, lower 
than the observed values of temperatures at sampling 
stations. First, tho simulated temperature values -are 
aveiage of temperature distribution in a day, while observed 
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Fig.5’8 Diurnal stream temperature 
simulation results 




River miles 
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values arc gral samples normally taken xn the afternoon. 
Doscrcpancy in those results is of the order of 5 to 10 °P. 
Secondly, local climatological data were collected in Baroda 
city end some of the sampling sites are 50 to 100 miles 
away from the metoreological station. 

Plow Augmentation: 

Dissolved soxygen concentrations WiTC falling “below 4.0 mg/l 
after river mile 70. Plow augmentation was carried out 
from Kadana reservoir to upgrade the water quality in Mahi 
river. It was found that minxmuin of 30 cft/sec. flow must 
bo discharged from Kadana reservoir. 

It was anticipated that after panam project on panam 
tributary is complete and com3.ssionod , the flow in panam 
river will reduce from 150 cft/sec. to 50 cft/sec. during 
summer months. Water quality degrades beyond limit in 
these conditions. These anticipated dissolved oxygen 
profiles using simulation model are not presented. Por this 
anticipated conditions, flow augmentation was explored. 

The flow augmentation results are presented in Table 4. 

Treatment Efficiency: 

As explained earlier, biochemical oxygen demand or bio- 
degradable organic matter content in the wastes from mint 
river donot have any significant effect on the dissolved 
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Table 4; Plow Augmoiitation Results. 


Plow rate from 
_cft ./sec . 

Kadana Panam 

Reservoir Reservoir 

Minimum 

DO c on-' 
cent rat ion 
mg/l 

location of 
minimum DO 
river mile 

Minimum 

conen. 

mg/l 

10.0 

50.0 

0.0 

11.0 

55.6 

28.1 

68.1 

1.4 

6.0 

27.7 

38.0 

78.0 

3.4 

5.0 

17.3 

42.65 

82.65 

3.8 

4.5 

14.8 

45.25 

85.25 

4.0 

4.5 

13.9 


oxygen concentrations of Mahi rivor. It is tiie high con- 
centrations of atnnonia nitrogen which consumes considerable 
amount of dissolved oxygen and DO sag is observed beyond 
mile. 70. It was found that 1,0 percent ammonia stripping 
from the mini river waste would nob allow the dissolved 
oxygen in river to fall below 4.0 mg/l.' 

Baroda Effluent Channel; 

As discussed earlier, this channel was proposed to carry 
waste from industries which are discharging, at present, 
their liquid waste into mini river, to the estuarian region 
of Mahi rivor, Big. 5.10 shows the anticipated dissolved 
oxygen profile after, this channel is commissioned. It was 
found that 30 percent of mini river waste can easily be ' 
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discharged into Mahi river even during critical months 
of summer. Thus, 30 percent reduction in designed capacity 
of Baroda effluent channel will result in considerable 
saving in its construction cost. 



ChapoGr 6 


COKCLUSIOHS 


Simula bio n studies for water quality in river Mahi indicates 
that classical Streeter Phelps' model for dissolved 
oxygen can not be applied where industrial wastes are 
dischaT*g.ed into river. If more reliable field survey 
data are available and properly used in such simulation 
models, a largo saving in both time and money can be 
achieved. 

It is concluded that minimum flow of 45.25 ft.^/sec. from 
Kadana reservoir and that of 85.25 eft. /sec. from Panam 
reservoir should be maintained to keep water q.uality in 
Mahi river above accepted standards. 

M 

1.0 percent treatment for ammonia removal from ; 0 iini river 
waste can be exercised to upgra.de the water quality in 
river Mahi. 

30.0 percent reduction in design capacity of Baroda 
Effluent Channel can be achieved without affecting the 
water quality in river Mahi, 
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[iecommoiidation for future Works 

Cn the present study, all the state variables, which 
lefine the e'.ltemativos, like flow augmentation and treat- 
nont efficiency were studied separately. By making use of 
jinivariate methods of optimization, the best and optimal 
combination of these state variables can be computed using 
3uch simulation models. 
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APPENDIX A 


SOLUTION TBCmiQUE 


Equation 4 can "be written for a control volume or element, 
in tho stream system as shown in Pig. 2.1 for nonuniform 
hydraulics as, 



where, 



[(AD. 


- (ADj ,]/Vi 




+ Sirifi-l 




V_ 


+ Qx. cx. 


+ 


' Sj_’ 


A A X = volume of control element, ft.^, 

) = moan cross •-sectional area 
of tho control volume ft.^ 


& O 

(ADy r") . 1 = Total longitudinal dispersion of the 

JU o X. 

constituent, ft . ^/sec .mg/l or temperature, 
ft^/sec.“°P, on tho inflow side of the ccaatrol 
volume, 


(ADj^ = Total longitudinal dispersion of the 

constituent, ft, ^/sec.mg/l or temperature, 
f t . ^/sec . -^p, on the outflow side of the control 
volume , 

Q. 1 = Bate of flow into the control volume, ft.^/sec, 



A~2 


c 


i-l 



G 

Qx. 

Gx. 



r 


Concentration of tlie constituent, mg/l, or 
temperature, in the inflowing water, 

Bate of flow out of the control volume, ft.^/sec. 

Goncentratioji of the constituent, mg/l, 
local inflows or withdrawals, ft.^/sec, 
Concentration of tho constituent, mg/l, 
or temperature, in Qx, and 

Sources or sinks of a nonconservative 
constituent, mg/l, or temperature, °F 


The foim of the differential ^equation used in the numerical 
solution of the transport equa bion is similar to Eq. A-1, 
except written in a slightly different form, 


aoi 


3T‘ = t(ADj^ “ (ADj^ 


3C. 


Q,.-l c.~l ~ Q. c + Qx- cx. 

-f - - -r, . nr ^ ^ ^ S . ^ 


(A. 2) 


The res\iLting finite difference form of Equation (A. 2) is 

^ ] Cl + 


At 


/ , ^ s cm; 

^i+1 “ 


(A. 3) 



^rhero, 


+ At Qz^ cz^ (a. 4 ) 

= -5- (A^ + = volume of 

olemoiit ±, (A. 5) 

an (3 all other terms are as previously defined. All of the 
values on the right-hand side of Eq,n. (A. 3) are Imown at time 
step and all of those of the left- hand side are unloiowns at 
time stop n+1. The coefficients on the left-hand side can 
be given as , 


a. 


= (AD^ ) 


At p ^ 
L^i-1 V^AZ “ "^-1-1 


(A. 6) 




^i 


1.0 + + ^-^^1^1 V^AZ ^i V^' 


(A. 7) 
(A. 8 ) 


M ethod _o_f ^ _Sj)J^t ion : 


Equation (A. 3) represents a tridiagonal sot of linear 
equations for the solution of C^' for all i's. 

This can bo represented in matriz form as, 






An efficient method that readily lends itself to a computer 
solution of the sot of Bq_imtion (A.9). The method of 
solution is as follows; 

(l) Divide through the first equation in (A.9) hy 
to obtain 


( 2 ) 


= ^1 (A. 10) 

where, 

^1 '^l^^l ^1 ” ^l^^l (A. 11) 

Combine Equation (A. 10) and the first equation in (A.9)} 
to oliminato ap and the result is,, 




+ Wo 0, 
2 5 


n+1 


(A. 12) 



whore 5 


2 = 




and G, 


Z2 — 8-2 
b2 “ 9-2 ^^ 


(A. 13 


OonbiiiG Equation (A. 12) and bho third equation in 
(A. 9) to climinato and tho 2?esult is, 


= G^ 


(A. 14 


irlioro 


'■'3 = ®3 


— a^ G2 

F^~arir" 

3 3 


(A. 15 


Eroceod through bhc equations, eliminating a^ and 
shoring bho values of II ^ £ind G^ given hy, 


^^i = b.'»a*nr.”7 » ^ = 2,3 ,..., h 

1 1 i“l 


(A. 16 


Z - a. G. . 

n 11 1*“1 _ Q '7' 

b-i = , 1 _ , n 


(A . 17 


The last equation is solved for C, by 


n p 

^h = '^h 


(A. 18 


Ti+l n-l-P Ti-l-l 

Solve for *^h back substitution: 

^ *^1+1’ ^ h-2,h-3,..., 1. 



